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We report BELLE measurements of the exclusive cross sections for the processes e+e" -> D M± D*^, e + e~ -> 
DD, e+e" -> D°D~n + , the first observation of ^(4415) -> DD* 2 (2460) decay and new state, Y(4660), using ISR. 
In addition, another cluster of events at around 4.05 GeV/c 2 is reported. 



1. Introduction 

Exclusive e + e~ hadronic cross sections to final 
states with charm meson pairs are of special in- 
terest since they provide information on the spec- 
trum of J PC = 1 charmonium states above the 
open-charm threshold. Parameters of these states 
obtained from fits to the inclusive cross section p] 
are poorly understood theoretically [2]. 

Initial-state radiation (ISR) is proved to be a 
powerful tool for measurement of the e + e~ exclu- 
sive hadronic cross sections at -y/s smaller than 
the initial e + e~ center-of-mass energy (-E c .m.) at 
£?-factories. ISR allows to obtain cross sections in 
a broad energy range while the high luminosity of 
the £?-factories compensates for the suppression 
associated with the emission of a hard photon. 

Here we report the first observation of the 
new charmonium-like state, F(4660) 3J, cluster- 
ing structure near 4.05 GeV/c 2 [4], the first mea- 
surement of the exclusive cross sections for the 
processes e+e" -> £>(*)±D*=f [5], e + e" DD 
[6], e + e~ — > D°D~tt + and the first observation of 
V>(4415) DD* (2460) decay [7]. The data sam- 
ple corresponds to a large integrated luminosity 
collected with the Belle detector [5] at the T(45) 
resonance and nearby continuum at the KEKB 
asymmetric-energy e + e~ collider [9]. 

2. Recoil mass technique 

There are two ways of ISR event reconstruc- 
tion: partial or full. 



In the full reconstruction method we select 
e + e~ — > DDjisr signal events by reconstruct- 
ing both the D and D mesons. In general, the 
7i sr is not required to be detected; its presence 
in the event is inferred from a peak at zero in 
the spectrum of the recoil mass against the DD 
system. The square of the recoil mass is defined 
as: 

M? ec (DD) = (E c . m , - E Dl5 ) 2 - p 2 D -, 

where E D jy and VjfD are ^ ne c - m - ener gy an d mo- 
mentum of the DD combination, respectively. 

To select e + e~ — > D^ + D*~ 7; sr signal events 
we use the partial reconstruction method that 
achieves high efficiency by requiring full recon- 
struction of only one of the mesons, the 7i sr , 
and the slow 7r~[ ow from the other D*~ [10] . In 
this case the spectrum of masses recoiling against 
the Z5(*) + 7i sr system 

M rec ( J D(*)+ 7isr ) = ^c.m.-^(.) + 7 ) 2 -^ ( . ) + 7 

peaks at the D*~ mass. Here E D (*)+ 1 andp D (.)+ 7 
are the cm. energy and momentum, respectively, 
of the D^*- )+ 7 i sr combination. This peak is ex- 
pected to be wide and asymmetric due to the pho- 
ton energy resolution function and higher-order 
corrections to ISR. To disentangle the contribu- 
tions from different final states and to suppress 
combinatorial backgrounds, we use the slow pion 
from the unreconstructed D*~ . The difference 
between the mass recoiling against D^ + y lsr and 
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£)W + 7r slow 7isr (recoil mass difference): 

AAf rcc = M rcc p(*)+7 isr )-M rcc (^(*)+7r s 7 ow7isr ) , 

has a narrow distribution around the nominal 
mD*- — ni-pa value, since the uncertainty in 7; sr 
momentum partially cancels out. 

3. Observation of the significant enhance- 
ment at 4.05 GeV. 

We identify e + e _ — * J / 'ipir + ir~ % SI process by 
peak in the distribution on the recoil mass against 
the J/-07r + 7r _ combination; J /if) is reconstructed 
in J/tp — > e + e~ and J/ip — > modes. 

Fig. [T] shows the ir + ir~J/ip invariant mass dis- 
tribution in the region of 3.8 — 5.5 GeV/c 2 . There 
is a clear enhancement at 4.25GeV/c 2 similar to 
that observed by the BaBar Collaboration [TTj . 
In addition, there is a clustering of events around 
4.05 GeV/c 2 that is significantly above the back- 
ground level. 

An unbinned maximum likelihood fit is applied 
to the 7r + 7r _ J/'0 mass spectrum in Fig. [T] a). 
Since there are two clusters of events in the 
mass distribution, we fit it with two coher- 
ent Breit-Wigner (BW) resonance functions as- 
suming there is no continuum production of 
e + e~ — > 7r + 7r~ J/ip. There are two solutions with 
equally good fit quality. The masses ((4008 ± 
40±28 4 ) MeV/c 2 and (4247 ± 12^) MeV/c 2 for 
the first and second states, respectively) and 
widths ((226 ± 44 ± 87) MeV/c 2 and (108 ± 19 ± 
10)MeV/c 2 ) of the resonances are the same for 
both solutions; the partial widths to e + e~ and the 
relative phase between them are different. The 
interference is constructive for one solution and 
destructive for the other. The statistical signifi- 
cance of the structure around 4.05 GeV/c 2 is es- 
timated to be 7.4(7 and is greater than 5<7 in all 
of the fitting scenarios that are considered. 

4. First observation of y(4660) state. 

Similar analysis is done for the ir + ir~ip(2S)ji sr 
final state. 

Fig-HJ b) shows the ir + ir~4>(2S) invariant mass 
for selected tp(2S) events, together with back- 
ground estimated from the scaled ip(2S) mass 
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Figure 1. Fit to the a): n + T:~'J/tp b): 
■k + tt~ t]j{2S) mass spectrum with two coherent 
resonances. The curves show the best fit and the 
contribution from each component. The dashed 
curves are for solution I, and the dotted curves 
for solution II. 



sidebands. Two distinct peaks are evident, one 
at 4.36 GeV/c 2 and another at 4.66 GeV/c 2 . 

An unbinned maximum likelihood fit that in- 
cludes two coherent P-wave Breit-Wigner (BW) 
functions and a constant, incoherent background 
is applied to the tt + tt~ 'ip(2S) mass spectrum 
in Fig. [TJ b). The fit results in two solutions 
with equally good fit quality, masses ((4361 ± 
9 ± 9) MeV/c 2 for the first state, (4361 ± 9 ± 
9) MeV/c 2 for the second state) and widths ((74± 
15±10)MeV/c 2 and (48 ± 15 ± 3) MeV/c 2 ). The 
interference is constructive for one solution and 
destructive for the other. A statistical signifi- 
cance of 5.8(7 is obtained for the peak around 
4.66 GeV/c 2 . 
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5. Measurement of the near-threshold 

e+e~ — > D^+D*^ cross section 

For the measurement of the exclusive cross sec- 
tion we determine the £)(*)+£)*" mass. The 
e+e- -> £>(*)+£)*- cross sections are extracted 
from the D^+D*~ mass distributions after back- 
ground subtraction using the relation described 
in [5]. The resulting exclusive e + e~ — » £>(*)+£>*- 
cross sections are shown in Fig. [5] The shape of 



to the BaBar data (blue points). The observed 
e + e~ — > £>Z3 exclusive cross sections are consis- 
tent with recent BaBar measurements [14] and 
are in qualitative agreement with the coupled- 
channel model predictions of Ref. [15]. This in- 
cludes a peak at 3.9GeV/c 2 that is seen both in 
Belle and BaBar mass spectra. 
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Figure 2. The exclusive cross sections for a) 
e+e" -> D*+D*- and b) e+e" -> D+D*-. 
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the e + e~ — > D* + D*~ cross section is complicated 
with several local maxima and minima. Aside 
from a prominent excess near the -0(4040), the 
e + e~ — > D + D*~ cross section is relatively fea- 
tureless. The measured cross sections are com- 
patible [T2] within errors with the D^D* ex- 
clusive cross section in the energy region up to 
4.260 GeV measured by CLEO-c [13]. 

6. Measurement of the near-threshold 
e + e~ — > DD cross section 

The e+e- -> D°i5 and e+e" -> D + D~ ex- 
clusive cross sections, measured with full event 
reconstruction method are shown in Fig. [3] Belle 
results, shown with the red points, are compared 



Figure 3. The exclusive cross sections for e + e — > 
D°D° and e+e" -> D+D~. The data are taken 
from the Durham database based on [6] and [14"] . 



7. Observation of 0(4415) -> DZ^(2460) de- 
cay 

We use the full reconstruction method de- 
scribed above to select e + e~ — * D a 'D~ tt + 7 j $ r 
signal candidates. The e + e~ — > D a D^n + 
cross section extracted from the background- 
subtracted D°D~ ~n + mass distribution demon- 
strates a prominent peak in a region of 0(4415) 
resonance. To study the resonant structure in 
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Figure 4. (a) The M£ > 0£>- 7r + spectrum for the 
Z?Z?2(2460) signal region. The threshold function 
is shown by the dashed curve, (b) The Mjjo D - n + 
spectrum outside the 13 (2460) signal region. 
The dashed curve shows the upper limit on the 
0(4415) yield at 90% C.L. Histograms show the 
normalized contributions from Mjjo and M D - 
sidebands. 



-0(4415) decays, we select D°D~ir + combina- 
tions from a ±100MeV/c 2 mass window around 
the nominal 0(4415) mass [16]. We perform 
a separate study of e + e~ — ► Z?£>2(2460) and 
e+e- ->■ ^(-D7r) non ^, (2460) . The M D o D -^+ spec- 
trum for the D£>2(2460) signal interval is shown 
in Fig. BJa). A clear peak corresponding to 
0(4415) -> ZXD^ (2460) decay is evident near the 
DI^ (2460) threshold. We perform a likelihood 
fit to M D o D -n+ distribution with the DD* 2 (2A60) 
signal parametrized by an s-wave RBW function. 
The significance for the signal is obtained to be 
~ lOcr. The obtained peak mass 171^144,15) = 
(4.411 ± 0.007(stat.)) GeV/c 2 and total width 
r tot = (77±20(stat.)) MeV/c 2 are in good agree- 
ment with the PDG [R)] values, the recent BES 
results [Ij and predictions of Ref . [2] . 

8. Conclusions 

In summary, we presented the first observation 
of the new charmonium-like state, F(4660), sig- 
nificant enhancement near 4.05GeV/c 2 , the first 
measurement of the exclusive cross sections for 
the processes e + e~ — * £)(*) ± £)*T ; e + e - _> jjTd^ 



e + e — > D°D ir + and the first observation of 
0(4415) -> I?!}?; (2460) decay. The obtained ISR 
results are in a good agreement with the recent 
CLEO-c and BaBar measurements of the exclu- 
sive cross sections. 
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